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CRPL Exponential Reference Atmosphere 


B. R. Bean and G. D. Thayer 


The background and development of an exponential model of atmospheric radio 
refractivity, the CRPL Exponential Reference Atmosphere, is outlined. A set of ray 
tracings for the model is presented in the form of tables of refraction variables for the 


complete range of observed values of surface refractive index. 
accuracy of the ray tracing and tabulation methods is made for these tables. 


A detailed analysis of the 
The variables 


are presented as numbers between 1 and 10 multiplied by the appropriate power of 10, thus 


maintaining a maximum number of significant figures. 


The tables may be used for the 


solution of practical refraction problems involving elevation angle errors, range errors, and 


similar quantities. 


1. Introduction 


Some exponential models of atmospheric radio 
refractivity have been introduced recently to 
replace the standard 4/3 earth treatment for 
applications imvolving long-range transmission 
paths, forward scatter predictions, or radar 
tracking at high altitudes [1, 2, 3, 4, 5].1 One of 
these models, developed at the Central Radio 
Propagation Laboratory and called the “CRPL 
Exponential Reference Atmosphere,” has’ been 


adopted for use by the National Bureau ot 
Standards in predictions of refraction phenomena. 

In this monograph, the background of the 
exponential reference atmosphere will be outlined 
and a set of radio ray tracings presented. These 
ray tracings may be used for the solution ot many 
problems involving refraction by means of the 
exponential reference atmosphere. 


2. Background 


In dealing with refraction problems it is con- 
venient to use the radio refractivity defined by: 


N=(n—1) X10, 


rather than the radio refractive index n. The 
basic predictor of refraction effects used in the 
exponential reference atmosphere is the surface 
value of refractivity at the transmitting point, 
N,. The reasons for this become clear upon 
examination of the normal, or average, N struc- 
ture of the atmosphere. It is found that the 
gradient of N which, with the usual assumption 
of horizontal homogeneity, is always with respect 
to height, is, on the average, a function of the 
surface value, N,. This is illustrated in figure 1 
where values of AN, the difference between JN, 
and the value of N at 1 km above the surface, are 
plotted against N,. The data represent 888 
six-yr monthly means from 45 U.S. weather 
Stations representing many climatically and geo- 
graphically diverse locations. As shown, AN 
‘1s approximately an exponential function of N,, 
_the relationship being expressed by: 


—AN=7.32 [exp {0.005577N;,}], (1) 
which is a least squares fit of In (—AN) to N,. 


! Figures in brackets indicate the literature references at the end of this 


monograph. 


It has also been found that the average decay 
of N with height is very nearly exponential. 
This can be seen from inspection of figure 2 
where averages and extremes of N structure are 
shown for medium high elevations on a semi- 
logarithmic plot. 

The exponential reference atmosphere was 
defined simply as that family of N profiles having 
a simple exponential decay with height and 
passing through the values NV, at the surface and 
N,+AN at a height of 1 km above the surface 
where AW is found from eq (1). N(h) is then 
given by: 

N=N, exp {—¢e(h—A;) }, (2) 
where h, is the surface elevation, h is the altitude 
above mean sea level, and c, is the decay constant 
given by: 


Ns 


aN (3) 


Ce=In 


Table 1 shows values of N,, AN, ¢., and h, 
for those profiles which were used in making up 
the tables of refraction variables given in ap- 
pendix B (a more complete table of this sort 
may be found in appendix A). The values of h, 
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Ficure 1. Siz-year mean AN as a function of 8-year mean N, at 0300 and 1500 G.m.t. for 45 U.S. weather stations. i 


were chosen to reflect the average elevation de- 
pendence of N,. 

It is found that several important advantages 
accrue from this particular choice of model, 
among which are: 

(1) The N profiles involved prove to be a fairly 
good representation of average N-structure over 
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Figure 2. Means and extremes of radio-refractivity height 
structure, and a comparison with extreme 4/3 earth 
profiles. 
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the first 3 km above the surface, a critical range i 
for refraction effects at low angles to the horizon. 

(2) The resulting relationship between the, 
N-gradient near the surface and JN, is found to be | 
very close to the observed average relationship | 
(for 5-yr mean profiles), being within +3 percent. 
over the normally encountered range of N, 
(from 300 to 365 for sea level, temperate zone | 
stations). This agreement is a requirement | 
for the successful prediction of average refraction | 
effects at elevation angles close to the horizontal. ‘ 

(3) The resulting profiles and their associated | 
gradients are continuous functions of height and | 
are in reasonably good agreement with the actual 


TaBLE 1. Exponential reference atmosphere profiles 
contained in appendix B 


} Effective 
Nz he AN Ce earth’s radius 
factor,* k 
i 
N ft N km-! 
200. 0 10,000 | —22. 3318 0. 118400 1.17779 
252. 9294 5, 000 —30. 6000 . 126255 1. 25551 
313.0 700 | —41. 9388 . 143859 1. 40229 
344. 5244 0 —50. 0000 _. 156805 1. 52479 
377. 2161 0 | —60. 0000 . 173233 1. 71320 
404. 8565 0 —70. 000U . 189829 1. 95921 
450. 0 0 —90. 0406 . 223256 2.77761 \ 
289. 0 0 | —36. 684830 0. 13574771 1. 3332410 
| 
i 1 dn ‘ 
*k=>————_ (notethat -—<0). To=a+h, where a is taken as 
To (dn ) dh | 
14+ ( — 
Ns dh h, 


6373.0150 km (3,960 miles). 


jobserved refractivity at all heights at frequencies 
‘sufficiently high so that the refractivity is not 
\influenced by the ionization in the ionosphere. 

_ An illustration of the effectiveness of the expo- 
nential reference atmosphere in predicting refrac- 
‘tion effects is given by figure 3, where the values of 
‘angular ray-bending are plotted against height for 
‘a ray starting at zero elevation angle. The agree- 
ment between the average bending for a large num- 
‘ber of profiles and the average of the values pre- 
‘dicted by the exponential reference atmosphere is 
‘seen to be quite good. Note however that the 
‘bending versus height predicted by the 4/3 earth 
model is markedly in disagreement with the 
observed for most heights. 

_ The preceding background material is intended 
only as an outline of the development of the expo- 
nential reference atmosphere. The remainder of 
‘this monograph will be devoted to an analysis of 
ray tracing theory and its application to the expo- 


nential reference atmosphere with particular 
emphasis on the preparation of the tables contained 
in appendix B. 
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FiGuRE 3. Average radio-ray bending for 22 5-yr. mean 
radio-refractivity profiles as a function of height, com- 
pared with the average of the values predicted by the exponen- 
tial reference atmosphere for each profile and the values 
predicted by the 4/3 earth treatment. 


3. Basic Theory 


_ The geometry of the radio-ray-refraction prob- 
‘lem is illustrated in figure 4, and the variables used 
in the following development are defined in the 
‘same figure. The path of the radio ray is the 
curved line R, which will also be used to represent 
the length of the ray path from the souree at 7) to 
‘the point at 7. If the transmitter is a radar set, 
the observed elevation angle to the target at 7 
would be 6, the true elevation angle would be 8B; 
thus the elevation error angle, ¢, is given by: 
{} 


e=H—B. (4) 


The angle 7 is a basic nieasure of refraction 
effects, and will be called the bending since it 
represents the angle through which the ray has 
_ been bent from its original direction due to refrac- 
tion. It can be seen from inspection of figure 4 
that: 
| 6—-hh=¢o—-7T, (5) 
from purely geometric considerations. 
_ Tf it is assumed that the refractive index is 
spherically stratified with respect to the earth, 
then Snell’s law (for the ray path) may be ex- 
_ pressed as: 


MP COS O=Np"p COS O9>= constant, (6) 
or 


| (+NX10-°)(a+h) cos 6 
| =(1+N,X10-5)(a-+hs) cos %, (7) 


where the radial distance from the center of the 
earth, r, equals a + h. The zero or s subscripts 
in (7) refer to the initial values at the earth’s 
surface. Equation (7) may be used to calculate 
the local elevation angle, 6, at any point along the 


ray and thus affords a complete description of 
the ray path. If: 


Nog COS O2=N, 1; COS 4), 


and 
M=N% + An, 
Tg=7,+Ar, 
and 
6.=6,+ A6, 


where An, Ar, and Aé@ are infinitesimals, then 
Mr, COS 6,=(n,+An) (7;+Ar)cos (6;+A8), 


and in the limits as Ar approaches zero this 
becomes: 


Mr, cos 6,=(n,+dn) (7, dr) (cos 6,—d6 sin 0,). (8) 


If (8) is expanded and products of differentials are 
omitted, then one obtains the differential equation: 


dn , dr 
tan Gd0 aay 


se cot pie cot 6. 
n r 


Noting that dr cot 6/r=rd¢/r=d¢, and, from 
geometry, di=d¢—dr, one obtains the classic ex- 
pression for the bending of a radio ray [6, 7]: 
dr= awe cot 6, 
n 


or 


n2 ad 
rua=— = cot 6. (9) 


Figure 4. The geometry of the radio-ray refraction problem. 


Note that (9) defines downward bending as a 
positive quantity. 

If the value of 7 has been obtained for some 
particular point (7, n) then the angle « may be 
found by the following scheme: 


Mr COS B=Np"o COS A 
from (6), and: 
r cos (6+7—€)=7p COS (A—Ee). 


(This is the equation of a straight line along the 
slant range, fo.) Combining these two equations, 
employing sum and difference trigonometric 
formulae, and solving for e will yield: 


cos r—sin 7 tan 6—n/n, 


e=Are tan - (10) 


n ‘ 
ae tan 6)>—sin 7—cos 7 tan 6 


Ss 


A similar expression for the angle 6 may be derived: 


n : 
——cos r—sin 7 tan 6, 
6=Are tan 


a0) 


: n 
sin r—cos 7 tan bot tan 0 


The distance along the earth’s surface from 7 to 7 | 


is found by: 
d=rop=1)(6—O% +7). 


The true distance from 7 to 7, called the slant 
range, Ro, is found from the law of sines or cosines: 


_rsin ¢ 


Ro= care Vr2+r?—2rirs COs ¢- 


(13) 


The range along the ray path, called the geo- 


metric range, Ff, is found by integrating over the. 


differential path lengths: 


R T 
r={ pe 
0 


, 
rg Sin 6 


while the effective or electrical (actually radio) 


range, R,., is found by: 


R T 
Re iD ndR= | i 
0 rg SIN 0 


The radio range error, AR,, and the geometric 
path length difference, AR, are defined as: 


AR, =R,— Ro, 


(15) 


and: 
AR=R—R.. 


It should be remembered that the above | 


equations are subject to the two usual restrictions 
of ray tracing: 


(1) The refractive index does not change ap- 


preciably in a wavelength. 

(2) The fractional change in the 
be small in a wavelength. 

Condition (1) will not be violated by the 
exponential reference atmospheres since there 
are no discontinuities in dn/dr.2 Condition (2) 
should be satisfied if: 


= < 0.002 fee» 


which would restrict the ray tracings for the ! 
exponential reference atmosphere to frequencies | 


over about 110 ke (for N,=450, for N,=200 the 
lower limit would be about 55 ke). 


2 This assumes also that | a | <i *, and 4 > 0. 


(12) | 
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(a4) | 


(16) 


(17) | 


spacing | 
between neighboring rays (initially parallel) must 
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| 4. Numerical Ray Tracing 


Because of the great amount of. calculation 
“necessary to obtain a set of ray tracings for the 
reference atmosphere it was advantageous to 
resort to digital computing machinery. A complex 
ray-tracing program had been written for the 
IBM 650, utilizing the floating point mode (8 
digits plus exponent), which was designed to 
handle arbitrary refractive index profiles having 
up to 100 layers and it was decided to analyze the 


exponential reference atmosphere with this pro- 


gram. This section will be devoted to an analysis 
of the ray tracing techniques employed in the 
machine program. 

The utilization of digital computer techniques 
for this type of ray tracing necessitates certain 
changes and alterations in the equations of 
refraction outlined in the previous section. 

To begin with, the basic accuracy of the ray 
tracing depends upon the accuracy with which the 
local elevation angles can be calculated using 
Snell’s law. In order to determine an angle in the 
neighborhood of 1 milliradian (mr) to an accuracy 
of 8 significant figures from the cosine of that 
angle, it is necessary to know the value of the 
cosine to 14 figures (cos r=0.999999xxxxxxxx, 1.€.: 
the 6 nines are nonsignificant) or to 1 part in 107‘. 

Rather than employ the relatively slow double- 
precision techniques required for the. machine 
evaluation of eq (6) an alternative formulation of 
Snell’s law was obtained by writing: 


Nolo 


— cos 
mr 


1—cos @=1— 


and taking advantage of the fact that 2 sin? 6/2=1— 
cos 6, the following equation was derived: 


a ; To 5 9 Gey 
6=2 Arc sin { [2 sin’ 5 sae 


Sa cos a (18) 


This equation is exact, and if each of tbe 
variables employed is correct to 8 figures, then 6 
will have an error no greater than 5 parts in the 
eighth figure. Since in performing ray tracings at 


initial elevation angles from 0 to 900 mr the value 


| 


of theta will never be as large as 1 radian, the 


Maximum error in calculating @ by this method 


should be 5 10~-* radians, or in general, not more 
than about 10~° percent of the elevation angle. 

With an accurate method of calculating the local 
elevation angles at various points along the ray 
path, the primary remaining problem is to com- 
pute accurately the total bending of the ray to 
each selected point. For this purpose a profile of 
N versus height is entered into the machine. 
The local elevation angle for each of the profile 
points may be calculated by means of eq (18), and 
if it is assumed that 


dN= constant 


dh (19) 


for each layer defined by succeeding pairs of 
points of the N profile, it can be shown that: 


Ben cot 61,91 (20) 
Ni, 2) 
where: 
O11 peta, 
j 2 


and the only approximation is made by letting: 


27, 2 2 2 
if an cot 6= i) dn cot @6== : f dn cot 6. 
1 7 M7N Ji M,2) J 1 


Since ;—n, never exceeds about 8010-* be- 
tween layer heights in the profiles, the maximum 
error entailed by this approximation cannot be 
very large (it is actually less than 10-*% of Ar in 
the worst case). 

It is evident from this method that 7 is essen- 
tially the weighted summation of the profile N 
gradients, where the cot 6; are the weighting fac- 
tors. It is important to note that, since for small 


1 see 
6 cot Oa the weighting factors are very large 


for the first few layers if 6=0. It is thus ex- 
pected that refraction effects at very-low eleva- 
tion angles would be extremely sensitive to the 
initial gradient of N. 

The integral for + was thus computed as the 
sum of Az for each layer: 


™ 
i—> Ati, (21) 
i=1 


where the only error of any magnitude is that 
made by assuming dn/dh=constant for each 
layer. In order to minimize this error it is neces- 
sary to restrict the thicknesses of each layer of 
the profile such that: A7;< constant, where the 
constant is some arbitrary angular amount. For 
the profiles presented in this paper the layer 
thicknesses were chosen such that A7 was always 
less than 2 mr for the first increment, and less 
than 3 mr for all succeeding increments. 

The approximate size of Ar, was found from the 
relation A7;+AN,/@, where 6, is the approximate 
value of 0, These were calculated for @=0, in 
which case 6,=~7j/Vkro/2 where k is the effective 
earth’s radius factor for the surface gradient and 


h is the mean height of the layer for which Az; is 
to be calculated. 
Thus: 


pve < constant. 


1 


It is thus possible to represent the given profile 
adequately with about 25 or less layers of constant 
N-gradient. it is also possible to calculate bend- 
ing using the approximation: n=constant for 
each layer. However it has been found that in 
order to insure sufficient accuracy about 1,000 
layers must be used to approximate the profile 
with this method [8]. 

The calculation of the distance is 
forward: 


straight 


di;=T0 (0;—40+ 71) (22) 


where 7 refers to the value of each variable at the 
top of the ith layer. 


(N,—N) 


6=Arc tan 


<10-*+[1— cos 7] —sin 7 tan 4 


vial 
ee ile 


Calculations of the elevation error angle are for) 
convenience made to depend on calculations of the 
angle 6: i 

e+7—6, 


since 6 is never larger than 1/2 7 and is thus more! 
easily calculated than e«. (Over the first few in-! 
crements of height the ray path is approximately) 
a circular arc, and thus e=6= 1/2 7, however as the’ 
range increases without limit the elevation angles’ 
formed by the ray path and the slant range ‘“‘path’’) 
both approach 7/2 and thus 60. This, of course, |! 
implies that lim e=hm e=T since r=e+6.) Equa-| 

oe my.) 


tion (11) for 6 is modified to: 


: Dye 
sin 7— cos 7 tan bo+7 tan 6 


for accuracy in calculation. 
The slant range is calculated with the law of 
cosines formulation (13), reduced to: 


Ro=¥ (r—10)?+2[1— cos 4] (25) 


again for accuracy of calculation. 

The geometric path length along the ray is 
calculated by the summation of the incremental 
path length through each layer (which will be 
denoted as As to avoid confusion with the range 
error, AR): 


R=>> As, (26) 
1=1 


which in turn is calculated by a law of cosines 
formulation similar to eq (25): 


As,=¥(ri—1i-1)?+2[1—cos (¢—¢;-))]. (27) 


Since the ray path through each layer is very 
nearly a circular arc, the maximum error in com- 
puting As; from (27) (assuming Az is exact) is: 


. =) 
2 sin 5 
error (As;) = E aes 100% 


2 
ee _ Ati 


ne e244 


X100% (to 7 sig. fig.), 


and since, as stated earlier, Av; is never as large 
as 0.003 radians, error (A,;) <0.0000375%, and 
hence: error (R) <0.00039%=2.5 m maximum 
since for the ‘‘worst case’ i.e.: N,=450, (the 
errors discussed in this and following sections 
are generally a maximum for the profile with 
the largest N gradient, in the case of the expo- 
nential reference atmosphere, the N,=450 profile), 
6=0, total r=31.5 mr, there would be about 10 


6 


increments of 3 mr each (if calculated with the} 
minimum number of layers), and Ry=650 km‘! 
when 7=31.4 mr, essentially the total amount) 
of bending. Actually there are some 30 layers’ 
calculated when 7=31.4 mr, and the calculated’) 
maximum error in # is only 0.0000514 percent: 
of 650 km or about 0.334 ms. It is felt that! 
calculation of ranges to an accuracy of +0.5 ms} 
or less is sufficient for any valid use to which the} 
exponential reference atmosphere might be put.) 
It should also be borne in mind that the accuracy) 
to which the distance along earth’s surface, the| 
slant range, the geometric path length, and the) 
radio range, can be calculated is also dependent) 
on the accuracy to which the angle ¢ has been 
calculated; this point will be discussed further) 
in the section dealing with the accuracy of the) 
tables. 1 

The radio range, (or “electrical path length’’,; 
“effective path length’’, etc.), is calculated simply | 
as:? | 


i 


™ ; i 
1 — > NG (28) \ 
i=1 | 


where:? 
asi= | ndR=7i,ds, UTM Ag, (29) | 
i=l 


Assuming that As; has been calculated accurately, 
the error in eq (29) for As’; depends on the extent 
to which: 


1 ; _ MAN -1 
RRs ee 
Since n;-;_; never exceeds 80 X10~°,and the error 
(defined as the difference between the above | 
terms) will never be as large as 2 percent of | 
N;N,-1, the maximum error in eq (29) should be 


less than 1.6X10~* percent. The worst case 
(again N,=450, @=0) yields a calculated maxi- 
“mum error of 2.810~* percent (at &,=800 km) 
or about 0.224 m. Note that the weighted 
average of n is always smaller than the arithmetic 


cm error in about 800 km, or 1.37 10-> percent). 
The errors discussed in the above paragraphs 

are, of course, much larger percentage-wise when 

applied to the values of the range errors: 


,mean of the end-point values, Ns and n,-;; this AR=R—R, 
indicates that the errors due to this approximation 
>will be of opposite sign to those due to using AR.=R.—Rp 


jeq (26) and (27) to calculate R. Since these 


,two errors will thus oppose each other when in the case previously mentioned as representing 


‘calculating F,, this quantity will tend to be more 
‘accurate than AR, the calculated error for the 
_N,=450, 6=0 case being only —0.110 m (11 


the largest errors (NV,=450, %=0, R)>800 km) 
the percentages are: Error (AR) + —0.67 percent 
maximum; error (AR,) = —0.06 percent maximum) 


5. Accuracy of the Tables 


In the preceding section on numerical ray- 
tracing the errors inherent in the equations and 
approximations employed were discussed at some 
length. It is important now to deal with those 
errors peculiar to the particular mode of computa- 
tion used in developing the tables, 1.e.: those 
errors which could be avoided by employing a 


-more accurate method of calculation using the 


same equations and approximations. The fol- 


: lowing material is a rather detailed analysis of 


: 


1 


j 


‘values used in the refractivity profiles. 


errors; for a concise presentation of the maximum 
errors likely to occur for each variable the reader 
is referred to the table at the beginning of ap- 
pendix B. : 

The basic source of error of this type is derived 
from the accuracy in the calculations of ou N 

he 
profiles for the exponential reference atmosphere 
were calculated to an accuracy of +5107" in 
n (or +5X10-° in N units). Thus the quantity 
N,—WNo.01, the difference between surface N and 
N at 10 m above the surface can be calculated to 


an accuracy of only about 5 significant figures 


(this is the worst case}, resulting in the following 
maximum errors: Error (@)<.+0.000047 mr 6,<4 
mr; error (@)<+0.000005 mr 6,>4 mr; and 


error in Ar<.0.0001 mr (from this cause only). 
_ These are random rounding errors and will not 


accumulate over a large number of height incre- 
ments. For 6, larger than 10 mr the accuracy 
will be unaffected by the mentioned errors in N, 
with the result that the largest error in 4 (for 4, 
of 400 and 900 mr) would be 3 or 4 units in 


error (Ar,)= (aro)| (1a Ge 


c. IN Ah; 


the eighth place: Error (@)<0.000004 mr 10 
mr<6@,<100 mr; error (@)<.0.00004 mr 6,>100 
mr. Hence the maximum error expected to be 
encountered in @ at any point is less than 510-8 
radians. 

The accuracy with which 7 can be calculated 
depends primarily on the accuracy with which 
N;,—WN;_; is known, so that 7 can be relied upon 
to 5 or 6 figures, depending on whether or not it 
is over 10 mr. Except for the first increment at 
6,=0, an error in the calculated value of Ar for 
some particular height increment due to a rounding 
error in N,—WN;_, will in general be at least % 
offset by an error of opposite sign over the fol- 
lowing height increment. An error in the first 
increment at @=0 will not usually be made up 
jater due to the excessively large weighting factor 


given the error by cot 6 in that case. Upon in- 
vestigating the average error in 7 at any height 
caused by the above mentioned errors, it is found 
that: partial error (7)<.0.0003 mr. 

This is an additional error in 7 due to the differ- 
ence between the true exponential decay of N and 
the linear approximation of N versus height as- 
sumed by the ray-tracing procedure. Errors of 
this type will be a maximum for the first increment 
of 7 for 6>=0, because of the rapid change of cot 4 
near 6=0. An equation has been derived for the 
difference in the incremental bending for the first 
layer at 6=O0 between the linear and exponential 
profiles (having the same value of N, and WN at 
7—?ro—Ah,), which is the error in Ar, due to the 
linear approximation: 

Wy 


CAh, 


sa! ce(Ahy)? _ 
3 10 : 


42 


+ 


error due solely to difference in N gradient assuming same cot @ 


(N,—N,)N, (c.— 


519249 O -59 -2 


N.—M 
N. 


© (Ewan) 


) 


[5-3] (Ah,)?? radians (30) 
BR 


error due solely to difference in cot @ where the 
subscript 1 refers to values of the variables at the 
first height increment, r—7)>=Ah,. In the case of 
the N,=450 profile the partial error in Av, from 
the first part of expression (30), amounts to only 
0.000760 mr, while the remainder of the error in 
Ar, from the second part of eq (30) amounts to 
only 0.00000725 wr, or only 1 percent of the first 
part. These errors drop rapidly with decreasing 
initial gradient of N, so that for the N,=377.22 
profile the error in Ary is only 0.000287 mr, or 
only about % as much as for the N,=450 profile. 

Errors in the value of the elevation error angle, 
e, will be of the same order of magnitude as those 
in 7 at any point, the method of computation being 
more precise than the values of 7 used in the 
calculation. 

Errors in the distance, d, and the ranges, Ro, RF, 
R, will be determined by the accuracy to which the 
angular distance can be calculated, 


$=0—8+7, 


which in turn depends on the accuracy to which 
6 and 7 are known. Since @ will have a maximum 
error of about 4107-5 mr and 7 a maximum error 
of about 8107‘ mr, the distance and ranges are 
calculated at least to an accuracy of about +5 m. 
However, if it is assumed that the profile is exact 
as represented by linear JN layers, the accuracy of 
the distance and ranges increases, with a maximum 
error of about +0.6 m. Actually the ranges for 
the N,=377.22 profile are accurate to about +1.8 
m, for the V,=200 profile the maximum error is of 
the order +0.6 m (about the same as if the profile 
were actually composed of linear N layers). 

The preceding estimates of maximum errors are 
based on the 4=0 ray; for any ray at 4 >0 the 
errors will in general be no larger than the limits 
imposed by computational accuracy (i.e.: errors 
not caused by the difference between ee ocnentil 
and linear N layers). Note that for very large 
6) the bending of the ray path is so slight that no 
significant differences can be detected between the 
arc lengths along the ray path, i.e.: the geometric 
range, and the slant range. In these cases the 
reported AR values are almost entirely one to four 
parts of the eighth significant figure of R, inter- 


f 
f 
| 


spersed with values given as exactly zero. These) 
are due to random rounding errors, and are often| 
negative. In general, to be considered as signifi-| 
cant (or real) a value of AR should be at least as. 
Jarge as two parts in the seventh figure of R. An) 
exception to this rule is for 6=900 mr; in this! 
case the angle ¢ can only be calculated to two | 
significant figures at the low heights, since 6—6@ 
is about 107° radians and @=6@ is about 1 radian, | 
and as a result of internal inconsistencies in the 
triangles formed by the ¢, and 7,, ry combinations, 
erroreous AF values as large as 0.002 percent of 
Ry are calculated. However, in none of the ray- 
tracings is there an entirely erroneous value of | 
AR or AR, reported which is larger than a few; 
centimeters (almost all such values are in the milli-; 
meter class, and differences between Ry and R as. 


\ 
; 


small as 107° km, or 0.001 mm are found in some}. 
of the ray-tracings for @=200 and 400 mr). A}, 


good rule to follow in using the values of AR and | 
AR, given in these ray-tracings would be to con-; 
sider that any of the values may be subject to an. 
error of at least +10 cm, and (in the case of the; 
larger values of AR) possibly as much as +50 cm. 

The values shown as “DELTA H” in the tables | 
are the corrections necessary to correct the height. 
given by a simple 4/3 earth atmosphere for the, 
distance given (assumed exact) to the height given | 
(also assumed exact). As a result, the values of; 
DELTA H have an error approximately equal to: 


error (delta h)= + (absolute erior in d) 5 — since the 


distances given are not really exact. The values 
shown as “DELTA THETA” are the corrections 
necessary to correct the elevation angle given by 
a simple 4/3 earth atmosphere at the given distance 
(assumed exact) to the value calculated for the 
exponential profile. This quantity is found by: 


ae=r—£ 


where 7 and ¢ are as calculated for the exponential 
reference atmosphere at the point where Aé is 
desired. The maximum error in this quantity 
therefore depends mostly on the errors in 7, hence: 
error (A@)<10~ radian. 


6. Tables of the CRPL Exponential Reference Atmospheres 


The tables themselves are simply a presentation 
of ray tracing variables. In these tables all 
angles are given in milliradians and all distances 
in kilometers. For each different profile (V,= 
200 through 450) ray tracings are given for 6 
from 0 to 900 mr, in approximately a geometric 
progression (i.e.: 0, 0.5, 1, 2, 4, 8, 15, 30, 65,100, 
200, 400, and 900 mr), at fixed ‘heights along the 
ray path. From left to right: first is given the 
SURFACE N, or N,, to 4 figures, not rounded. 


Exact values of N, used are given in table 1. 
The next column is the value of 6, or “INIT (ial) 
THETA”, in milliradians. The “HEIGHT” 
is given in kilometers. 
are given with fixed decimal points, the remaining 
quantities are listed as floating point numbers in 
the following fashion: +X.XXXX +E where the 
first sign refers to the sign of the number, 


X.XXX which is to be multiplied by 10 to the | 


power +£ depending on the second sign. 


The preceding quantities | 


! 


i 


| 


} 


' THETA at same distance, 


Examples: 


0.0000 = zero 

1.2345 = PBs Cx<ilder 

4 eo — 122345110? 
—1.2345 3=—1.2345 x10? 

1.2345—-3= 1.234510? 


—1.2345—3=— 1.2345 107°. 


The following quantities are listed as floating 
point numbers in this manner: “DELTA H”= 


' correction to 4/3 height at same distance to equal 
| (HRIGHT””’, 


in kilometers; ‘‘THETA’’=@6 or 
local elevation angle, in milliradians; ‘“DELTA 
THETA” =correction to 4/3 earth @ to equal 
in milliradians; 


‘ “DISTANCE”=7¢ or distance along earth’s 
' surface, in kilometers; ‘““TAU’’=7, angular bend- 


ing, in milliradians; “ERROR ANGLE’’=.e, 
elevation error angle, in milliradians; “SLANT 


| RANGE” —&,, straightline distance from ray 


origin to point at fA, d, in kilometers; “DELTA 
R”=AR, difference between distance from h=0 


to given height along ray path and slant range, in 


kilometers; “DELTA R—E’’=AR., difference be- 
tween radiopath length and slant range, in 
kilometers. 

The following ray-tracing variables have not 


: been printed, but can be easily calculated from 
_ those given, as follows: 


pee. ~) EDISTANCE” 
1% “INITIAL RADIUS”’ 


where 7>=6373.0150-+-A, (kilometers) 
g=7—e=“TAU”—“ERROR ANGLE”, 


R=R,+AR=“SLANT RANGE” 
+“DELTA R”, 


R,=R,+AR,=“SLANT RANGE” 
+“DELTA R—E”. 


Values of elevation of the ray starting-point, h,, 
are given in table 1. 

The last profile, for N,=289.0, is a special 
profile which was prepared for the purpose of 
making direct comparisons between 4/3 earth 
refraction and an exponential profile having the 
same WN gradient at the surface. It is intended 
for applications where it is desired to replace 
the 4/3 earth treatment with a single standard 
atmosphere having the advantage of differing 
very little from the 4/3 earth near the surface. 


The authors wish to acknowledge the technical 
assistance of Miss B. A. Cahoon, and Mrs. B. J. 
Weddle, and particularly of K. A. Norton whose 
help in working out some of the mathematical 
methods employed was invaluable. 


7. References 


fi] B. R. Bean and G. D. Thayer, On models of the 
atmospheric radio refractive index,- Proc. IRE 47, 
No. 5, 740 (May 1959). 

(2] B. R. Bean, Atmospheric bending of radio waves, 
eee presented at International Conference on 

lectromagnetic Wave Propagation, Liege, Belgium, 
October, 1958 (to be published by the Academy 
Press, spring 1959). 

[8] L. J. Anderson, Tropospheric bending of radio waves, 
Trans. Am. Geophys. Union 39, No. 2, 208 (February 
1958). 

[4] J. R. Bauer, W. C. Mason, F. A. Wilson, Radio refrac- 
tion in a cool exponential atmosphere, Lincoln 


Laboratories, Technical Report No. 186 (August 
27, 1958). 

[5] G. Millington, Propagation at great heights in the 
atmosphere, The Marconi Rev. 21, No. 131, 143, 4th 
Quarter (1958). 

[6] M. Schulkin, Average radio-ray refraction in the lower 
atmosphere, Proc. IRE 40, 554 (May 1952). 

[7] W. M. Smart, Spherical astronomy, Chapter III 
(Cambridge University Press, London, 1931). 

[8] G. H. Millman, Tropospheric effects on radar target 
measurements, Proceedings of the Seventh Weather 
Radar Conference, held at the Deauville Hotel, 
Miami Beach, Florida (November 17-20, 1958). 


8. Appendix A. Tables of Profile Variables for the CRPL Exponential Reference 


Atmosphere 


The following equations were used in determin- 
ing the various profile variables: 


—AN=7.32 exp {0.005577 N,} 


fo Ne 
Come INGEN 


—dN,=c.N, 


1 


== Ns=1+ N,X10-§8, 
lay (¢.N,) X10>° 


k 


ry. is taken as 6373.014987 km for all values of 
N, in these tables. 

The table giving N,, —AN, c., and —dN, as 
functions of k was prepared using an iterative 
method for solution of the transcendental equa- 
tions involved. In all these tables the accuracy 
may be taken as +1 in the last digit listed. 


TaB_eE A-1l. 

N, a —A N Co = dN; 0 k 

200 | 22.3317700 | 0. 118399435 | 23. 6798870 1. 17769275 
210 | 23.6125966 | .119280212 | 25. 0488444 1. 18991401 
220 | 24.9668845 | . 120458179 | 26. 5007993 1. 20315637 
230 | 26.3988468 | .121916361 | 28. 0407631 1. 21752719 
240 | 27.9129385 | . 123642065 | 29. 6740955 i. 23314913 
250 | 29.5138701 | . 125626129 | 31. 4065323 1. 25016295 
260 | 31. 2066224 | . 127862319 | 33. 2442030 1. 26873080 
270 | 32.9964614 | . 130346887 | 35. 1936594 1. 28904048 
280 | 34.8889558}  . 133078254 | 37. 2619112 1. 31131073 
290 | 36. 8899932 | . 136056720 | 39. 4564487 1. 33579768 
300 | 39.0057990 |. 139284287 | 41. 7852861 1. 36280330 
310 | 41.2429556 | . 142764507 | 44. 2569972 1. 39268608 
320 | 43. 6084233 |. 146502381 46. 8807620 1. 42587494 
330 | 46. 1095611 . 150504269 | 49. 6664087 1. 46288731 
340 | 48. 7541501 . 154777865 | 52. 6244741 1. 50435338 
350 | 51.5504184 | . 159332141 | 55. 7662495 1. 55104840 
360 | 54.5070651 |. 164177379 | 59. 1038565 1. 60393724 
370 | 57. 6332884 | . 169325150 | 62. 6503054 1. 66423593 
380 | 60. 9388149 | 174788368 | 66. 4195799 1. 73349938 
390 | 64. 4339281 | . 180581312 | 70. 4267116 1. 81374807 
400 | 68.1295015 | . 186719722 | 74. 6878887 1. 90765687 
410 | 72. 0370324 | . 193220834 | 79. 2205420 2. 01884302 
420 | 76. 1686780 | . 200103517 | $4. 0434770 2. 15232187 
430 | 80. 5372922 | . 207388355 | 89. 1769927 2. 31525447 
440 | 85.1564647 | . 215097782 | 94. 6430240 2. 51823286 
450 |  90.0405683 | . 223256247 | 100. 4653113 2. 77761532 
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TABLE A-2 
—AN N, Co —dNo k 
20 | 180. 226277 | 0.117626108 | 21.1993155 | 1.15617524 
22 | 197.316142 | .118216356 | 23.3259953 | 1.17457412 
24 | 212.917967 | .119594076 | 25.4637276 | 1. 19366808 
26 | 227.270255 | .121491305 | 27.6113599 | 1.21348565 
28 | 240. 558398 | .123746115 | 29. 7681671 | 1. 23406110 
30 | 252.929362 | .126255291 | 31.9336703 | 1. 25543336 
32 | 264.501627 | .128950180 | 34.1075325 | 1.27764560 
34 | 275.372099 | .131783550 | 36. 2895127 | 1.30074523 
36 | 285.621054 | .134721962 | 38. 4794288 | 1.32478398 
38 | 295.315731 | .137741207 | 40.6771452 | 1.34981825 
40 | 304.513148 | .140823306 | 42. 8825481 1. 37590934 
42 | 313.261483 | .143955014 | 45. 0955611 1. 40312414 
44 | 321.602888 | .147125889 | 47.3161107 | 1. 43153539 
46 | 329.573439 | .150328075 | 49.5441405 | 1.46122250 
48 | 337.204713 | .153555418 | 51.7796106 | 1. 49227226 
50 | 344.524418 | .156803056 | 54.0224815 | 1.52477960 
52 | 351.557000 | .160067149 | 56. 2727266 | 1. 55884863 
54 | 358,324138 | .163344614 | 58.5303179 | 1.59459364 
56 | 364.845143 | .166633002 | 60.7952415 | 1.63214058 
58 | 371.137293 | .169930326 | 63.0674811 | 1. 67162830 
60 | 377.216108 | .173234984 | 65.3470266 | 1.71321044 
62 | 383. 095581 . 176545680 | 67. 6338699 | 1. 75705732 
64 | 388.788373 | .179861358 | 69.9280046 | 1. 80335830 
66 | 394.305974 | .183181171 | 72.2294300 | 1. 85232456 
68 | 399.658845 | .186504431 | 74.5381454 | 1.90419925 
70 | 404.856538 | .189830583 | 76.8541525 | 1. 95922635 
72 | 409.907798 | .193159183 | 79.1774555 | 2.01772514 
74 | 414.820650 | .196489873 | 81.5080567 | 2. 08002556 
76 | 419.602477 | .199829385 | 83.8459677 | 2. 14650999 
78 | 424.960086 | .203156494 | 86.1611914 | 2. 21761358 
80 |} 428.799768 | .206492043 | 88.5437400 | 2. 29383429 
82 | 433.227348 | .209828917 | 90.9036251 | 2. 37574437 
84 | 437.548229 | .213167031 | 93.2708570 | 2. 46400458 
86 | 441.767432 | .216506335 | 95.6454475 | 2. 55938222 
88 | 445. 889634 | .219846812 | 98.0274147 | 2. 66277367 
90 | 449.919193 | .223188453 | 100. 4167688 | 2. 77523207 
92 | 453.860184 | . 226531281 | 102. 8135290 | 2. 89800399 
94 | 457.716416 | . 229875327 | 105.2177108 | 3.03257531 
96 | 461.491458 | .233220637 | 107. 6293319 | 3. 18073184 
98 | 465.188659 | .236567271 | 110.0484114 | 3.34463902 
100 | 468.811163 | .239915290 | 112.4749663 | 3. 52694820 
TaBLe A-3 
k N, —AN Ce —dNo 
1.0 0.0 0.0 0.0 0.0 
1.2 | 217.689023 | 24.6471681 | .120160519 | 26. 1576259 
1.3 | 275.037959 | 33.9367000 | .131692114 | 36. 2203304 
1.4 | 312.297111 | 41.7747176 | .143600133 | 44. 8459068 
1.5 | 339.003316 | 48.4839018 | .154339490 | 52. 3215987 
1.6 | 359.298283 | 54.2941700 | .163827653 | 58. 8629945 
1.7 | 375.341242 | 59.3759008 | .172203063 | 64. 6349115 
1.8 | 388.391792 | 63.8586055 | .179626805 | 69. 7655765 
1.9 | 399.243407 | 67.8426334 | . 186242834 | 74. 3562234 
2.0 | 408.424907 | 71.4070090 | . 192172034 | 78.4878451 
2.1 | 416.304322 | 746148487 | .197514185 | 82. 2260091 
2.2 | 423.146728 | 77.5171828 | .202351472 | 85. 6243634 
2.3 | 429.148472 | 80.1557288 | .206751820 | 88. 7272277 
2.4 | 434.458411 | 82.5649192 | .210771674 | 91.5715264 
2.5 | 439.191718 | 84.7734613 | .214458304 | 94. 1883108 
2.6 | 443.438906 | 86. 8054237 | .217851443 | 96. 6038056 
2.7 | 447.272272 | 88.6811886 | . 220984823 | 98. 8403838 
2.8 | 450.750273 | 90.4181120 | .223887193 | 100. 9172131 
© 523. 299600 | 135.5109159 | .299693586 | 156. 829534 
4/3 | 289.036274 | 36. 6922523 | .135758874 | 39. 2392391 
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9. Appendix B. Tables of Refraction Variables for the Exponential Reference 
Atmosphere 


The following table of estimated maximum errors should serve as a guide to the accuracy of the 


tables. 
Errors in elevation angle, 0: 


6.<4 mr +0. 00005 mr 


4 mr<6)<100 mr + .000005 mr }nearly independent of N,. 


6) >100 mr + .00004 mr 


Errors in 7, ¢ (in milliradians): 


N.= 450 404. 8 377. 2 344. 5 313 252. 9 200 
4=0 +0. 001 0. 00065 0. 0005 0. 0004 0. 0003 0. 0002 0. 0002 
=1° . 0003 . 00015 . 0001 . 00008 . 00006 . 00005 . 00005 
4=3° . 00004 . 000025 . 00002 . 000017 . 000015 . 000013 . 000012 
Errors in Ry, R, R., or Ah (in meters): 

ve 450 404. 8 377. 2 344. 5 313 252. 9 200 
=0 +5. 0 72, Uf 18 1.2 0.8 0. 65 0. 6 

A > 1° 0. 4 0. 3 0. 25 0. 2 ING . 15 . 14 


Assume that error in AR or AR, is 0.5 percent or 


+0.1 m, whichever is larger. 
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U. S. GOVERNMENT PRINTING OFFICE : 


THE NATIONAL BUREAU OF STANDARDS 


The scope of activities of the National Bureau of Standards at its major laboratories, in Washington, 
D.C., and Boulder, Colo., is suggested in the following listing of the divisions and sections engaged 
in technical work. In general, each section carries out specialized research, development, and engi- 
neering in the field indicated by its title. A brief description of the activities, and of the resultant 
publications, appears on the inside front cover. 


WASHINGTON, D.C. 


Electricity and Electronics. Resistance and Reactance. Electron Devices. Electrical Instru- 
ments. Magnetic Measurements. Dielectrics. Engineering Electronics. Electronic Instru- 
mentation. Electrochemistry. 

Optics and Metrology. Photometry and Colorimetry. Optical Instruments. Photographic 
Technology. Length. Engineering Metrology. 

Heat. Temperature Physics. Thermodynamics. Cryogenic Physics. Rheology. Molecular 
Kinetics. Free Radicals Research. 

Atomic and Radiation Physics. Spectroscopy. Radiometry. Mass Spectrometry. Solid State 
Physics. Electron Physics. Atomic Physics. Neutron Physics. Radiation Theory. Radio- 
activity. X-ray. High Emergy Radiation. Nucleonic Instrumentation. Radiological Equip- 
ment. 

Chemistry. Organic Coatings. Surface Chemistry. Organic Chemistry. Analytical Chemistry. 
Inorganic Chemistry. Electrodeposition. Molecular Structure and Properties of Gases. 
Physical Chemistry. Thermochemistry. Spectrochemistry. Pure Substances. 

Mechanics. Sound. Mechanical Instruments. Fluid Mechanics. Engineering Mechanics. Mass 
and Seale. Capacity, Density, and Fluid Meters. Combustion Controls. 

Organic and Fibrous Materials. Rubber. Textiles. Paper. Leather. Testing and Speci- 
fications. Polymer Structure. Plastics. Dental Research. 

Metallurgy. Thermal Metallurgy. Chemical Metallurgy. Mechanical Metallurgy. Corrosion. 
Metal Physics. 

Mineral Products. Engineering Ceramics. Glass. Refractories. Enameled Metals. Consti- 
tution and Microstructure. 

Building Technology. Structural Engineering. Fire Protection. Air Conditioning. Heating 
and Refrigeration. Floor, Roof, and Wall Coverings. Codes and Safety Standards. Heat 
Transfer. Concreting Materials. 

Applied Mathematics. Numerical Analysis. Computation. Statistical Engineering. Mathe- 
matical Physics. 

Data Processing Systems. SKAC Engineering Group. Components and Techniques. Digital 
Circuitry. Digital Systems. Analog Systems. Applications Engineering. 

« Office of Basic Instrumentation. ¢ Office of Weights and Measures. 


BOULDER, COLORADO 


Cryogenic Engineering. Cryogenic Equipment. Cryogenic Processes. Properties of Materials. 
Gas Liquefaction. 

Radio Propagation Physics. Upper Atmosphere Research. Jonosphere Research. Regular 
Prediction Services. Sun-Earth Relationships. VHF Research. Radio Warning Services. 
Airglow and Aurora. Radio Astronomy and Arctic Propagation. 

Radio Propagation Engineering. Data Reduction Instrumentation. Radio Noise. Tropo- 
spheric Measurements. Trophospheric Analysis. Propagation-Terrain Effects. Radio-Mete- 
orology. Lower Atmosphere Physics. 

Radio Standards. High-Frequency Electrical Standards. Radio Broadcast Service. Radio and 
Microwave Materials. Electronic Calibration Center. Microwave Physics. Microwave Circuit 
Standards. 

Radio Communication and Systems. Low Frequency and Very Low Frequency Research.. 
High Frequency and Very High Frequency Research. Modulation Systems. Antenna Research. 
Navigation Systems. Systems Analysis Field Operations. 
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